Afterlive: A performant code for Vlasov-Hybrid simulations

A parallelized implementation of the Vlasov-Hybrid method [Nunn, 1993] is presented. This method is a hybrid between a gridded Eulerian description and Lagrangian meta-particles. Unlike the Particle-in-Cell method [Dawson, 1983] which simply adds up the contribution of meta-particles, this method does a reconstruction of the distribution function $f$ in every time step for each species. This interpolation method combines meta-particles with different weights in such a way that particles with large weight do not drown out particles that represent small contributions to the phase space density. These core properties allow the use of a much larger range of macro factors and can thus represent a much larger dynamic range in phase space density. The reconstructed phase space density $f$ is used to calculate momenta of the distribution function such as the charge density $\rho$. The charge density $\rho$ is also used as input into a spectral solver that calculates the self-consistent electrostatic field which is used to update the particles for the next time-step. Afterlive (A Fourier-based Tool in the Electrostatic limit for the Rapid Low-noise Integration of the Vlasov Equation) is fully parallelized using MPI and writes output using parallel HDF5. The input to the simulation is read from a JSON description that sets the initial particle distributions as well as domain size and discretization constraints. The implementation presented here is intentionally limited to one spatial dimension and resolves one or three dimensions in velocity space. Additional spatial dimensions can be added in a straight forward way, but make runs computationally even more costly.Comment: Accepted for publication in Computer Physics Communication

PICPANTHER: A simple, concise implementation of the relativistic moment implicit Particle-in-Cell method

A three-dimensional, parallelized implementation of the electromagnetic relativistic moment implicit particle-in-cell method in Cartesian geometry (Noguchi et. al., 2007) is presented. Particular care was taken to keep the C++11 codebase simple, concise, and approachable. GMRES is used as a field solver and during the Newton-Krylov iteration of the particle pusher. Drifting Maxwellian problem setups are available while more complex simulations can be implemented easily. Several test runs are described and the code's numerical and computational performance is examined. Weak scaling on the SuperMUC system is discussed and found suitable for large-scale production runs.Comment: 29 pages, 8 figure

Physics case for an LHCb Upgrade II - Opportunities in flavour physics, and beyond, in the HL-LHC era

The LHCb Upgrade II will fully exploit the flavour-physics opportunities of the HL-LHC, and study additional physics topics that take advantage of the forward acceptance of the LHCb spectrometer. The LHCb Upgrade I will begin operation in 2020. Consolidation will occur, and modest enhancements of the Upgrade I detector will be installed, in Long Shutdown 3 of the LHC (2025) and these are discussed here. The main Upgrade II detector will be installed in long shutdown 4 of the LHC (2030) and will build on the strengths of the current LHCb experiment and the Upgrade I. It will operate at a luminosity up to 2×1034 cm−2s−1, ten times that of the Upgrade I detector. New detector components will improve the intrinsic performance of the experiment in certain key areas. An Expression Of Interest proposing Upgrade II was submitted in February 2017. The physics case for the Upgrade II is presented here in more depth. CP-violating phases will be measured with precisions unattainable at any other envisaged facility. The experiment will probe b → sl+l−and b → dl+l− transitions in both muon and electron decays in modes not accessible at Upgrade I. Minimal flavour violation will be tested with a precision measurement of the ratio of B(B0 → μ+μ−)/B(Bs → μ+μ−). Probing charm CP violation at the 10−5 level may result in its long sought discovery. Major advances in hadron spectroscopy will be possible, which will be powerful probes of low energy QCD. Upgrade II potentially will have the highest sensitivity of all the LHC experiments on the Higgs to charm-quark couplings. Generically, the new physics mass scale probed, for fixed couplings, will almost double compared with the pre-HL-LHC era; this extended reach for flavour physics is similar to that which would be achieved by the HE-LHC proposal for the energy frontier

LHCb upgrade software and computing : technical design report

This document reports the Research and Development activities that are carried out in the software and computing domains in view of the upgrade of the LHCb experiment. The implementation of a full software trigger implies major changes in the core software framework, in the event data model, and in the reconstruction algorithms. The increase of the data volumes for both real and simulated datasets requires a corresponding scaling of the distributed computing infrastructure. An implementation plan in both domains is presented, together with a risk assessment analysis

Determining Pitch-Angle Diffusion Coefficients for Electrons in Whistler Turbulence

Transport of energetic electrons in the heliosphere is governed by resonant interaction with plasma waves, for electrons with sub-GeV kinetic energies specifically with dispersive modes in the whistler regime. In this paper, particle-in-cell simulations of kinetic turbulence with test-particle electrons are performed. The pitch-angle diffusion coefficients of these test particles are analyzed and compared to an analytical model for left-handed and right-handed polarized wave modes

Analytical treatment of particle motion in circularly polarized slab-mode wave fields

Wave–particle interaction is a key process in particle diffusion in collisionless plasmas. We look into the interaction of single plasma waves with individual particles and discuss under which circumstances this is a chaotic process, leading to diffusion. We derive the equations of motion for a particle in the fields of a magnetostatic, circularly polarized, monochromatic wave and show that no chaotic particle motion can arise under such circumstances. A novel and exact analytic solution for the equations is presented. Additional plasma waves lead to a breakdown of the analytic solution and chaotic particle trajectories become possible. We demonstrate this effect by considering a linearly polarized, monochromatic wave, which can be seen as the superposition of two circularly polarized waves. Test particle simulations are provided to illustrate and expand our analytical consideration